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ABSTRACT 

Gravitational lcnsing provides a unique and powerful probe of the mass distribu- 
tions of distant galaxies. Four-image lens systems with fold and cusp configurations 
have two or three bright images near a critical point. Within the framework of sin- 
gularity theory, we derive analytic relations that are satisfied for a light source that 
lies a small but finite distance from the astroid caustic of a four-image lens. Using a 
perturbative expansion of the image positions, we show that the time delay between 
the close pair of images in a fold lens scales with the cube of the image separation, with 
a constant of proportionality that depends on a particular third derivative of the lens 
potential. We also apply our formalism to cusp lenses, where we develop perturbative 
expressions for the image positions, magnifications and time delays of the images in a 
cusp triplet. Some of these results were derived previously for a source asymptotically 
close to a cusp point, but using a simplified form of the lens equation whose validity 
may be in doubt for sources that lie at astrophysically relevant distances from the 



caustic. Along with the work of |Keeton et al. (20051, this paper demonstrates that 



perturbation theory plays an important role in theoretical lensing studies. 
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1 INTRODUCTION 



Gravitational lensing, or the bending of light by gravity, offers an exciting synergy between mathematics and astrophysics. 
Singularity theory provides a powerful way to describe lensing near critical points ( |Schneider et al.p~9 92 Pott ers et al.|200Tj |, 
which turns out to have important implications for astrophysics and the quest to understand dark matter. Four-image lensed 
quasars can be broadly classified, according to image geometry, as either folds, cusps or crosses. We are interested in folds 
and cusps, which occur when the light source is close to the caustic curve, along which the lensing magnification is infinite 
(see top row of Fig. [I]). Fold lenses feature a close pair of bright images, while cusp lenses have a close triplet of bright images 
(see bottom row of Fig. jlj. By expanding the gravitational potential of the lens galaxy in a Taylor series, one finds that the 
image magnifications satisfy simple analytic relations, viz. 

Ma|-|mb|~0 and \(ia\ - \^b\ + \ftc\ ~ (1) 

for fold pairs and cusp triplets, respectively ( Blandford fc Narayan|1986{|Ma o 1992, Schneider & Weiss 1992; Schneider et al. 
1992||Zahkarov|1995|rGaudi fc Petters|2002a|b Keeton et al. pool ' |2005 \ For a fold pair, the two images have opposite parity, 
hence the negative sign. For a cusp triplet, the middle image (B) has opposite parity from the outer images (A and C). Note 
that in practice, one works with the image fluxes, which are directly observable, rather than the magnifications, which are 
not. This leads to the equivalent relations: 

Fa-Fb „ , „ Fa-F b +Fc 



Riold = 



and -Rcubd = 



0, 



Fa + Fb ' P Fa + Fb+Fc 

where the image flux F, is related to the source flux Fq by Ft 



(2) 



\m\F . 
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Observationally, the fold and cusp relations are violated in several lens systems (e.g., Hogg & B landfordj 19 94 Falco et al. 



1997 Keeton et al.|il997 2003 20051. These so-called "flux-ratio anomalies" are taken as strong evidence that lens galaxies 



contain significant small-scale structure (e.g., Mao & Schneider 1998 Metcalf & Madau 2001 Chiba 2002). Since global 



modifications of the lens potential (modeled with multipole terms, for example) cannot explain the anomalous flux ratios 
( Evans fc Witt|2003 Congdon fe Keeton|2005 Yoo et al.;2005, 20061, it is generally believed that lens galaxies must contain 
substructure in the form of mass clumps on a scale > 10 6 Mq ( Dobler fe Keeton||2006[). Indeed, substru cture is observed in 
at least two lens systems, in the form of a dwarf galaxy near the main lens galaxy |Schechter fc Moore 1993. McKean et al. 



20071. The substructure could well be invisible, though; numerical simulations of structure formation in the cold dark matter 



(CDM) paradigm predict that galaxy dark matter halos are filled with hundreds of dark subhalos (e.g., Klypin et al. 1 1999 



Moore et al. 19991. The abundance of flux-ratio anomalies in lensed radio sources implies that lens galaxies contain ~2% 



(0.6-7% at 90% confidence) of their mass in substructure ( |Dalal fc Koch anck 2002]) . This may be somewhat higher than the 
amount of CDM substructure (e.g., Mao et al.|2004 1, although new higher-resolution simulations are attempting to refine the 
theoretical predictions and make more direct comparisons with lensing observations (e.g., Gao et al.|2004 Strigari et al.|2007| 



Diemand et al.|[2007 l. While there is still work to be done along these lines, it is clear that violations of the ideal fold and 



cusp relations provide important constraints on the small-scale distribution of matter in distant lens galaxies. 



Keeton et al. ( 2003 2005 I pointed out that the ideal fold and cusp relations only hold for a source asymptotically close to 



the caustic. If we want to use flux-ratio anomalies to study small-scale structure in lens galaxies, it is vital that we understand 
how much J?f id and i? C usp can deviate from zero for realistic smooth lenses just because the source lies a small but finite 
distance from the caustic. To that end, Keeton et al. ( 2003 1 used a Taylor series approach to demonstrate that 7? C us P = to 
lowest order, and used Monte Carlo simulations to suggest that R C us P oc d 2 , where d is the distance between the two most 
widely separated cusp images. In order to derive the leading non-vanishing term in 7? C usp analytically, it would be necessary to 
extend the Taylor series to a higher order of approximation. This has not been done before, since including higher-order terms 
substantially complicates the analysis. Instead, previous authors have made assumptions about which terms are important 
and which terms can be neglected, in order to obtain an analytically manageable problem. 

As we shall see, perturbation theory provides a natural way to overcome the difficulties of the Taylor series approach. 



Keeton et al. ( 2005 I used perturbation theory to show that for fold lenses, i?foid is proportional to the image separation d\ of 



the fold pair. We extend their analysis in several important ways. We derive the leading-order non-vanishing expression for 
^cusp using perturbation theory. We also show, for the first time, how the combination of singularity theory and perturbation 
theory can be used to study time delays between lensed images for both fold and cusp systems. Working in analogy to flux 
ratios, we derive time-delay relations for fold and cusp lenses, which we will apply (in a forthcoming paper) to observed lenses 
in order to identify "time-delay anomalies." As time delays join flux ratios in probing small-scale structure in lens galaxies 
(see Keeton & Moustakas 2008 I , singularity theory and perturbation theory once again provide the rigorous mathematical 
foundation. 



2 MATHEMATICAL PRELIMINARIES 



To study lensing of a source near a caustic, it is convenient to work in coordinates centered at a point on the caustic. 
Nonspherical lenses typically have two caustics. The "radial" caustic separates regions in the source plane for which one 
(outside) and two (inside) images are produced^ Within the radial caustic is the "tangential" caustic or astroid, which 
separates regions in the source plane for which two (between the two caustics) and four (inside the astroid) images are 
produced. We are interested in four-image lenses, where the source is within the astroid. We therefore make no further 
reference to the radial caustic. A typical astroid is shown in Figure [T] for a lens galaxy modeled by a singular isothermal 
ellipsoid (SIE), which is commonly used in the literature (e.g., Kormann et al.|[l994 \ and appears to be quite a good model 
for real lens galaxies (e.g., Rusin fc Kochanek|2005 Koopmans et al.|2006[ ). It is customary to define source-plane coordinates 
(2/1,3/2) centered on the caustic and aligned with its symmetry axes. However, for fold and cusp configurations, where the 
source is a small distance from the caustic, it is more natural to work in coordinates (ui, 112) centered on the fold or cusp point. 
For convenience, we define the u\ axis tangent to (and the 112 axis orthogonal to) the caustic at that point. Transforming 
from the (3/1,3/2) plane to the (111,112) plane requires a translation plus a rotation. To derive this coordinate transformation, 



we follow the discussion in Appendix Al in Keeton et al. (2005 I, which summarizes the results of Petters et al. (2001 1 



We begin by considering the lens equation y = x-V^(x), which maps the image plane to the source plane. The solutions 
to this equation give the image positions x = (xi,X2) corresponding to a given source position y = (3/1,3/2)- The function 
ip(x) is the scaled gravitational potential of the lens galaxy projected onto the lens plane. A caustic is a curve along which 
the magnification is infinite, i.e., det(<9y/<9x) = = 0, where <9y/<9x is the Jacobian of y, and is known as the inverse 



1 We are neglecting faint images predicted to form near the centers of lens galaxies, because they are highly demagnified and difficult to 
detect. 



Analytic Relations for Fold and Cusp Lenses 3 



magnification matrix. We choose coordinates such that the origin of the source plane (y = 0) is on the caustic. In addition, 
we require that the origin of the lens plane (x = 0) maps to the origin of the source plane. We are interested in sources that 
lie near the caustic point (y = 0), which give rise to lensed images near the critical point (x = 0). In this case, we may expand 
the lens potential in a Taylor series about the point x = 0. We then find that the inverse magnification matrix at x = is 
given by 



<9x 
where 



1 - 2a 
-b 



1 - 2c 



a = -V>n(0), 



6 = ^12(0), 



c = 2^22(0). 



(3) 



(4) 



The subscripts indicate partial derivatives of ip with respect to x. Note that ip has no linear part (since y = when x = 0). 
For y = to be a ca ustic point, we must have (1 — 2a)(l — 2c) — b 2 = 0. In addition, at least one of (1 — 2a), (1 — 2c), and 
must be non-zero (Petters et al. 



2001 



p. 349). Consequently, (1 — 2a) and (1 — 2c) cannot both vanish. Without loss of 



generality, we assume that 1 — 2a 7^ 0. 

We now introduce the orthogonal matrix (see Petters et al.|2001 p. 344) 

1 I" l-2o -S 

M = 

y/(l-2a) 2 + b 2 I b 1 - 2a 
which diagonahzes <9y/9x| . We then define new orthogonal coordinates by 

= (61,62) = Mx, u= (ui,u 2 ) = My. 



(5) 



(6) 



Note that the coordinate changes are the same in the lens and source planes. The advantage of using the same transformation 
in both the lens and source planes is that the lens equation takes the simple form 



u = 6 Vip(0), 

and that the inverse magnification can be written as 



det 



(-) 

\do) 



(7) 



(8) 



The old and new coordinate systems in the source and image planes are shown in Figure [T] Since the caustic in the source 
plane maps to the critical curve in the image plane, the origin of the (81,62) frame is determined from that of the (111,112) 
frame. The orientation of the (61,82) axes is determined by the matrix M, and is not necessarily related to the tangent to 
the critical curve. 

Using the local orthogonal coordinates u and 6, Petters et al. (2001 p. 346) showed that x = is a fold critical point if 
and only if the following conditions hold 

(l-2a)(l-2c) = S 2 , l-2a#0, ^222(0) ± 0. (9) 
For a cusp, the third condition above is replaced by the requirements that 

^222(0) = 0, Viaa(O) / 0, V>222 2 (0) / 0. (10) 

Note in particular that 7/^222 (0) = for a cusp while 1/^222 (0) 7^ for a fold; this indicates that these two cases must be treated 
separately. 

We are interested in obtaining the positions, magnifications and time delays of images near critical points. Since these 
quantities depend only on the behavior of the lens potential near a fold or cusp point, we can expand ip(0) in a Taylor series 
about the point 6 — 0. To obtain all the quantities of interest to leading order, we must expand the lens potential to fourth 
order in 6 ( |Petters et al.|2001| pp. 346-347): 



^(81,62)= 2(1 



K) 



?i + - 6 2 2 +edl + f6l$2+gdi8l + h8l + k8 i 1 +m( 



'l$2 + n8l$l +p6i( 



i + r8t, 



(11) 



where the coefficients {K, e, f, g, . . . ,r} are partial derivatives of the potential evaluated at the origin. Lensing observables 
are independent of a constant term in the potential, so we have not included one. Since 6 = maps to u = 0, any linear 
terms in the potential must vanish. In the second order terms, the coefficients of the 8182 and 8 2 terms are set to and 1/2, 
respectively, in order to ensure that the point = is a critical point (see Appendix Al of Keeton et al.|2005). 



3 THE FOLD CASE 



In this section, we use perturbation theory (e.g., Bellman 19661 to derive an analytic relation between the time delays in 
a fold pair. To derive this expression, we must first obtain the image positions at which the time delay is evaluated. These 
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Figure 1. The curves show the caustic (top) and critical curve (bottom) for a singular isothermal ellipsoid (SIE) lens with minor-to- 
major axis ratio q = 0.5. The local orthogonal coordinates defined by the rotation matrix M are indicated for a fold point (left) and 
a cusp point (right). In the top two panels, a source (triangle) with position (1/1,3/2) measured from the center of the caustic (astroid) 
has position (111,112) in the rotated coordinates centered on a fold point (left panel) or a cusp (right panel). In the bottom two panels, 
the four images (triangles) have positions (x\,X2) measured from the center of the critical curve (ellipse) and positions (81,82) in the 
rotated coordinates. 



results were derived by Keeton et al. (2005 I. We offer a summary of their analysis in Section 3.1 and present our new results 



for the time delay in Section [5\2 



3.1 Image Positions 

Since we are considering a source near a fold point, we write its position in terms of a scalar parameter e, which we take to 
be small but finite. In particular, let u — » eu. Combining equations |7| and ( |11[ ) we can write the lens equation as 



eiti 



£U2 



KOi- (3e0? + 2/0i0 2 + 



gel) 



(4k 9\ + 3m 6(62 + 2n Ml + 



VV2) 



(/ B\ + 2g M2 + 3ft el) - (m ef + 2n 0i 2 + 3p 0i02 + 4r &. 



(12) 
(13) 



Analytic Relations for Fold and Cusp Lenses 5 



(see Petters et al. 2001 Theorem 9.1). To find the image positions, we expand 0i and 02 in a power series in e. Since the 
left-hand sides of equations ( 12 I and (13 1 are accurate to 0(e), the right-hand sides must be accurate to the same order. 
Noting that the lowest-order terms on the right-hand side are linear or quadratic in 6, we write 

0i = ai e 1/2 + (3 1 e + (e) 3/2 . 
ot2 e 



02 = a 2 e 1/2 + fc e + O (e) 3/2 . 



(14) 
(15) 



Substituting into the lens equation, we obtain 

= (a 1 K)e 1/2 -(3aie + 2a 1 a2f + a 2 2 g-/3 1 K + u 1 )e + 0(e) :i/2 , (16) 
= — (a 2 / + 2ai«2.g + Sa^h + U2)e — [2ai/?i/ + 2(ai/?2 + cx2pi)g + 6«2/32ft + a?m + 2Q 2 Q2n + 3aia 2 p + 4afr] e 3 ^ 2 

+0(e) 2 . (17) 



Note that these equations are carried to different orders in e, since the leading-order term in equation ( 12 1 is linear in 6, while 
the leading-order term in equation ( |13[ ) is quadratic in G. 

Since e is non-zero, equations and ( 171 must be satisfied at each order in e. We can work term by term to solve for 
the coefficients ai, a.2 and 02, and then write the image positions as 

Sghm - g 2 u 2 _ , ^3/2 



3ftif 



Of 



-U2 1/2 

3ft 6 



9h 2 K 



e + 0(e) 



(18) 
(19) 



where the ± labels indicate the parities of the images. From these equations, we see that two images form near the point 
6 — on the critical curve, provided that (— M2/3/1) > 0. Since h ^ for standard lens potentials (e.g., an isothermal ellipsoid 
or isothermal sphere with shear), we must have U2 > 0. In other words, the source must lie inside the caustic in order to 
produce a pair of fold images. In practice, a more useful quantity is the image separation, given by 



di 



-4w 2 
3ft 



e V2 + Q ( e ) 



3/2 



(20) 



3.2 Time Delays 

To find the time delay between the two fold images, we begin with the general expression for the scaled time delay (e.g. 



Schneider et al.| 1992 1 



f(«) = :^ = lie-ui 2 -^). (21) 

To 2 
The scale factor is given by 

1 + z L D L Ds foo , 

where Dl, Ds and Dls are the angular-diameter distances from the observer to lens, observer to source, and lens to source, 
respectively. The lens redshift is denoted by zl- Making the substitution u — > (eui, €112), we have for the two fold images 



nO , - J_-i- (eU2 )3/2 + 0(e) 2 ; (23) 



t + =t{0+) = -J-—(eU2r^ + 0(ey. (24) 



The time delay between images is then (cf. Schneider et al.||l992 pp 190 - 191) 

(25) 



Af fold = f- - f+ = (eu2f 2 + O (e) 2 =-^dl + 0(, 



which is positive, in agreement with the general result that images with negative parity trail those with positive parity. We 
find that the only coefficient from the lens potential that enters the expression for the differential time delay is the parameter 
ft = ^222(0)/6. We also see that to leading order in e, the image separation and the differential time delay depend only on the 
«2 component of the source position. Unlike the image positions, our expression for the time delay does not involve any of the 
fourth-order terms in the potential. This is because the time delay involves the potential directly, while the image positions 
depend on first derivatives of the potential. This means that all fourth-order terms in the potential enter the time delay at 
O (e) 2 , while these same terms enter at O (e) in quantities involving derivatives. 
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To summarize, Aff id oc (ev,2) 3 ^ 2 oc df. For comparison, Rtoid oc di. Since d\ is small, a violation of the ideal relation 
Affoid = is more likely to indicate the presence of small-scale structure in the lens galaxy than would be indicated by a 
non-zero value of 7?f id. 

Our analytic expression for Affoid is only valid for sources sufficiently close to the caustic that higher-order terms are 
negligible. To quantify this statement, we compare our analytic approximation with the differential time delay computed 
numerically from the exact form of the lens equation. The numerical analysis requires a specific lens model, so we consider 
an SIE with axis ratio q = 0.5 as a representative example. We use the software of Keeton (20011 to compute the astroid 
caustic, and then choose a point on the caustic, far from a cusp, to serve as the origin of the (ui, U2) frame. For a given value 
of U2, we solve the exact lens equation to obtain the image positions and time delay for the fold doublet. The top left panel 
of Figure [2] shows Affoid in units of 6% as a function of U2 in units of Be, where 9e is the Einstein angle of the lens. The 
analytic and numerical results are in good agreement for sources within 0.05#e of the caustic, although the curves do begin to 
diverge as U2 increases. The difference between the numerical and analytic curves (which represents the error in the analytic 
approximation, denoted by e) is shown in the middle left panel. The bottom left panel shows the logarithmic slope of the 
error curve, d(lne)/d(lnu2). Together, the middle and bottom left panels verify that our analytic approximation is accurate at 
order e 3//2 . Furthermore, these panels suggest that the next non- vanishing term is of order e 5 ^ 2 and has a positive coefficient. 
The interesting implication is that the coefficient of the e 2 term seems to vanish. At our current order of approximation we 
are not able to determine whether this is rigorously true, and if so, how general it is; we merely offer the remark in the hope 
that it may be useful for future analytic studies. For now, we focus on the verification that our analytic approximation is 
accurate at the order to which we work. 

Since the source position 112 is not observable, we compare the analytic and numerical time delays as a function of the 
image separation d\ in the right-hand panels of Figure [2] The range of d\ corresponds to that used for 112 in the left-hand panels. 
For a canonical fold lens with di = 0A69e ( |Keeton et al.|[2005[ ), our analytic expression gives a very good approximation, 
indicating that our analysis can be applied in astrophysically relevant situations — not just when the source is asymptotically 
close to the caustic, but even when it lies some small but finite distance away. 



4 THE CUSP CASE 

We now apply our perturbative method to the case of a source near a cusp point. This approach has not been applied to 
cusp lenses before. Appendix A of Keeton et al. ( 2003 1 derives the image positions and magnifications for a cusp triplet 
assuming a simplified form of the lens equation. As we noted in Section [l] this simplified lens equation assumes that certain 
terms may be set to zero, using criteria that are less than rigorous. We use the lens equation derived from the fourth-order 
lens potential, and use perturbation theory to verify the results of Keeton et al. ( 2003 1 and extend the analysis to a higher 
order of approximation. We also study time delays for a cusp lens for the first time. Our analysis does not involve simplifying 
assumptions, and indicates that perturbation theory is a powerful method in the study of lensing. 



4.1 Image Positions 



We again expand the image positions, magnifications and time delays in the parameter e, but with one notable difference. 
For a cusp oriented in the u\ direction, a small "horizontal" displacement of etti from the cusp point permits a "vertical" 
displacement of only e 3 ^ 2 W2 (see Fig.Jl]), since larger vertical displacements would imply a source position outside the caustic 
( Blandford fc Narayan|1986 |. The lens equation is then 



em = K Oi - (3e6» 2 + 2/6>i t 



■gOl)- UkOf + 3m9 2 9 2 + 2n9 1 9l+p9l) 



3/2 

e ' U2 



(26) 

(/(9 2 + 2g9 1 9 2 ) - (m9'i + 2n9 2 1 92 + 3p9 1 9l+4r9^) (27) 

(see Petters et al.||2001 Theorem 9.1), where the 9\ term of equation (13 1 does not appear, since 1^222(0) = for a cusp, 
corresponding to h — in equation (111. 

As before, we write the image positions as a series expansion in e, but now keeping an additional term (i.e., 7;e 3//2 ). This 
is necessary since the vertical component of the source position enters the lens equation as e 3//2 «2, rather than 012 as in the 
fold case. We have 



1/2 

9 2 = Q 2 e 1/2 + (3 2 e + 72 e 12 + O (e) 
The lens equation then becomes 



= aiKe 1 ' 2 — (3ea 2 + 2fct\oi2 + gee 2 — (3\K + u 1) e 

- [Akal + 3ma 2 Q 2 + pa 2 + 2 ai (na 2 2 + 3ej3i + ffo) + 2a 2 (//3i + gfo) - Kji] e 3/2 + O (e) 2 



(28) 
(29) 



(30) 
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Figure 2. The top panels show the time delay for a fold pair as a function of source position (left) and image separation (right), 
for an SIE lens with axis ratio q = 0.5. The solid line shows our analytic approximation while the dotted line shows the exact result 
obtained by solving the lens equation numerically. The quantities 112 and di are defined in the text, and are shown here in units 
of the Einstein angle, 0g. The time delay is given in units of 0^. The middle panels show the error in the analytic approximation, 
e = Affoid (numerical) — Aff id (analytic), due to our neglect of higher-order terms. The bottom panels show the logarithmic slope of 
the error curve, rf(lne)/d(lnU2) on th e left and d(lne)/d(ln di) on the right. There is some numerical noise in the logarithmic slope due 
to numerical differentiation. Together the middle and bottom panels verify that our analytic expression is accurate at order e 3 / 2 (cf. 
equation |25[i . 



= —(/<*! + 2gQlQ2)e - 



/ 3 2 2 3 \ 

Ilt2 + ma x + 2na 1 ci2 + 3paia 2 + 4ra 2 + 2/aipi + 2ga 2 f5\ + 2ga\f32) 



3/2 



- {piiffc + 2gp 2 ) + al(3m/3i + 2nf3 2 ) + 3aa(p/3i + 4r/3 2 ) + 2(?a 27 i + 2ai [a 2 (2n/3i + 3p/3 2 ) + fji + 972]} < 
+0 (e) 5/2 . 



(31) 



As in the fold case (see equation 16 1, we find u\ = 0. Note that 72 appears only in the e coefficient of equation (31 1, but in 
a term multiplied by ai . Hence, it will not be possible to solve for 72 ■ Fortunately, it turns out that the expressions we would 
like to derive do not involve this parameter. 
We can now write the lens equation as 







(gal 



PiK + m) 1 



\P&2 



+ 2a 2 (/A + gP 2 ) - K lx ] e 3/2 + O {ef , 



(32) 
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■ 4ra 3 2 + 2ga 2 (3 1 ) e 3/2 - [Pl(f@i + 2g(5 2 ) + 3a 2 (p/3i + 4r/3 2 ) + 2 fl a 2 7i] <? + O (e) 5/2 . 



= - (U2 

Recalling that these equations must be satisfied at each order in e, we then solve for the unknown coefficients: 



0i 



7i = 



ff«2 + «1 

K ' 

(5/g 2 + 5gKp)a 4 2 + (6/gm + 3Kpui)a 2 2 + K 
2if (3g 2 a| + GKra\ + gm ) 
(/<? 3 - 2 3 2 iCp + 12/ffJTr + GK 2 pr)a s 2 + (2fg 2 Ul - 2gKpu 1 + 12fKr Ul )a 3 2 + fgu 2 a 2 



K 2 (3g 2 a 2 + QKrai, + gui) 



where a 2 satisfies the cubic equation 
gu\ t Ku 2 



a 3 2 + 



7«2 + 



0. 



2Kr + g 2 4Kr + 2g 2 

This is equivalent to equation (A8) from |Keeton et ah] ( |2003| ), after making the replacements oi 2 
—a/4. To leading order, the image positions can be written as 



z, K -> c, g 



ga 2 + ui 
= Pie = e, 



K 



(h 



a 2 e 



1/2 



(33) 

(34) 
(35) 
(36) 

(37) 
-6/2, r 

(38) 
(39) 



which are equivalent to equation (A7) from Keeton et al. (20031. The distance between any two images i and j is then 



di 



„(») 



+C(e) 



(40) 



4.2 Magnifications 

Substituting our perturbative expressions for the image positions into equation djjj, we find that the inverse magnification of 
a cusp image is given by: 

A* 1 = ~2[gui + 3{g 2 + 2Kr)al\ £ + [(24/r - 12flp)ajj + (6/Cp - ifg)oaPi ~ {8g 2 + 24Kr)a 2 f3 2 - gKji] e 3/2 + O (e) 2 . (41) 
We then find from equation |2]| that 



Rc 



+ O (e) « Tlcuspd , 



(42) 



where d = maxij dij is the largest separation between any two of the three cusp images. This expression shows that correction 
terms to the ideal cusp relation enter at second order in the image separation, which agrees with the conjecture of [Keeton] 
(20031. To see this, we define rat = l/i^ 1 ], which allows us to write 



et al. 



-Rcusp — 



m B mc — niAmc + mAm B 



(43) 



m B mc + rriAmc + m^ms 

If the leading-order term in the numerator vanishes, so does the leading-order term in 7? cusp . The zeroth-order term in i? CU sp 
corresponds to a term of O (e) 2 in the numerator, since the leading-order term in the denominator is O (e) 2 . By substituting 
the solutions for a 2 into the numerator, we find that i? CU sp = to lowest order, in agreement with Keeton et al. (20031. We 
repeat this procedure for the next-leading term of O (e) 5//2 in the num erator, and find that R C us V = at linear order in d [i.e. 
O(e) 1 ^ 2 ] as well; this result was unattainable using the formalism of 



Keeton et al. 



( 2003b. To proceed to higher order, we 



must extend our perturbative analysis by including terms of the form 6te 2 in equations (281 and (29 1 for the image positions. 
We denote the coefficient of d 2 by A CUBp , which we do not here write down, since that would require several pages. Given the 
complexity of this term, it is not practical to evaluate ^4 C us P analytically. However, we have numerically computed -Rcusp for 
the case of an SIE model with q — 0.5 (see Fig.[3| and find that R CUBp oc u\ oc d 2 . While it is conceivable that y4 cusp might be 
zero for some specific lens model, clearly this is not the case in general. We have thus demonstrated analytically that i? CUS p 



vanishes through linear order in d, placing the numerical result of Keeton et al. ( 2003 1 on solid mathematical ground 



4.3 Time Delays 

For the cusp case, the scaled time delay takes the form 



2 

We find 



•eui) +' 



3/2 



u 2 y] - 4:(ei, e 2 ). 



X - [(3g 2 + 6Kr)a 4 2 + 2gu x a\ + (K — l)u\] e 2 + O (e] 



5/2 



2K 



(44) 



(45) 
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Figure 3. -Rcusp as a function of u\ (left) and d (right) for an SIE with q = 0.5, obtained by solving the lens equation numerically. 



corresponding to a differential time delay of 



Af, 



(ij) 



1 

Ik 



23 (4* 



) mi + 3Ku 2 



(46) 



Unlike the fold case, the time delay for a pair of cusp images depends on both source coordinates (ui, W2). This means that it 
is not possible to write our current expression strictly in terms of observables, such as the image separation. Instead, all we 
can say is that the time delay scales quadratically with e, or alternatively, with the fourth power of the image separation. 

In the fold case, we found that the time delay scales as e 3 ^ 2 and only depends on the lens potential through the parameter 
h. For a cusp, however, h = 0, so it is not surprising that the lowest-order term in the time delay is of O (e) 2 . Furthermore, if 
we had not included the ji e 3//2 terms in our expansions of the image positions for a cusp [equations ( 28 



and 



291] 



it would not 

2 



have been possible to obtain a perturbative expression for the time delay; instead, we would simply have found f = + O (e) 



5 SUMMARY 

We have developed a unified, rigorous framework for studying lensing near fold and cusp critical points, which can (in principle) 
be extended to arbitrary order. We have found that the differential time delay of a fold pair assumes a particularly simple 
form, depending only on the image separation and the Taylor coefficient h — ^222 (0)/6. This result is astrophysically relevant, 
since it is quite accurate even for sources that are not asymptotically close to the caustic. We have also obtained perturbative 
expressions for the image positions, magnifications and time delays of a cusp triplet. These results rest on the key insight that 
a source at a given distance e from a cusp along the relevant symmetry axis of the caustic can only move a perpendicular 
distance of e 3 ^ 2 in order to remain inside the caustic (Blandford fc Narayan|l986 1. We have shown rigorously that the distance 



Keeton et al. 



( 2003) is correct. We have also demonstrated that 



dependence of the magnification ratio R C us P conjectured by 
the leading-order expression for the image positions is given by the relations presented by Keeton et al. (20031, and have 



provided the necessary framework for deriving the image positions corresponding to a Taylor expansion of the lens potential 
at arbitrary truncation order. Finally, we have derived cusp time delays analytically for the first time. Our results provide a 
rigorous foundation for identifying anomalous flux ratios and time delays in gravitational lens systems, and for using them to 
study small-scale structure in the mass distributions of distant galaxies. 



ACKNOWLEDGMENTS 



We thank A. O. Petters and Peter Schneider for helpful discussions and suggestions, and for their careful reading of the 
manuscript. We also thank the anonymous referee for useful comments. 



10 A. B. Congdon, C. R. Keeton and C. E. Nordgren 



REFERENCES 



Bellman, R. 1966, Perturbation Techniques in Mathematics, Engineering, and Physics (Mineola, NY: Dover) 
Blandford, R., & Narayan, R. 1986, ApJ, 310, 568 
Chiba, M. 2002, ApJ, 565, 17, arXiv:astro-ph/0109499 

Congdon, A. B., & Keeton, C. R. 2005, MNRAS, 364, 1459, |arXiv:astro-ph/0510232| 



Dalai, N., & Kochanek, C. S. 2002, ApJ, 572, 25, |arXiv:astro-ph/01 11456 



Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 657, 262 

Dobler, G., & Keeton, C. R. 2006, MNRAS, 365, 1243, |arXiv:astro-ph/0502436| 

Evans, N. W., & Witt, H. J. 2003, MNRAS, 345, 1351, |arXiv:astro-pIi/0212oT3| 

Falco, E. E., Lehar, J., & Shapiro, I. I. 1997, AJ, 113, 540, |arXiv:astro-ph/9612182| 

Gao, L., White, S. D. M., Jenkins, A., Stoehr, F., & Springel, V. 2004, MNRAS, 355, 819 



Gaudi, B. S., & Petters, A. O. 2002a, ApJ, 574, 970, arXiv:astro-ph/0112531 



— . 2002b, ApJ, 580, 468, arXiv:astro-ph/0206162 
Hogg, D. W., & Blandford, R. D. 1994, MNRAS, 268, 889, |arXiv:astro-ph/9311077l 
Keeton, C. R. 2001, |astro-ph/0102340l 

Keeton, C. R., Gaudi, B. S., & Petters, A. O. 2003, ApJ, 598, 138, |arXiv:astro-ph/0210318| 

— . 2005, ApJ, 635, 35, |arXiv:astro-ph/0503452l 
Keeton, C. R., Kochanek, C. S., & Seljak, U. 1997, ApJ, 482, 604, |arXiv:astro-ph/9610T63l 
Keeton, C. R., & Moustakas, L. A. 2008, arXiv:0805.0309, 805 



Klypin, A., Kravtsov, A. V., Valenzuela, O., & Prada, F. 1999, ApJ, 522, 82, |arXiv:astro-ph/9901240 

Koopmans, L. V. E., Treu, T., Bolton, A. S., Buries, S., & Moustakas, L. A. 2006, ApJ, 649, 599, |arXiv:astro-ph/0601628| 

Kormann, R., Schneider, P., & Bartelmann, M. 1994, A&A, 284, 285 

Mao, S. 1992, ApJ, 389, 63 

Mao, S., Jing, Y., Ostriker, J. P., & Weller, J. 2004, ApJ, 604, L5 

Mao, S., & Schneider, P. 1998, MNRAS, 295, 587, |arXiv:astro-ph/9707187| 



McKean, J. P. et al. 2007, MNRAS, 378, 109 
Metcalf, R. B., & Madau, P. 2001, ApJ, 563, 9 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, J., & Tozzi, P. 1999, ApJ, 524, L19, |arXiv:astro-ph/9907411| 
Petters, A. O., Levine, H., & Wambsganss, J. 2001, Singularity Theory and Gravitational Tensing (Boston: Birkhauser) 
Rusin, D., & Kochanek, C. S. 2005, ApJ, 623, 666, |arXiv:astro-ph/0412001| 
Schechter, P. L., & Moore, C. B. 1993, AJ, 105, 1 

Schneider, P., Ehlers, J., & Falco, E. E. 1992, Gravitational Tenses (Berlin: Springer- Verlag) 
Schneider, P., & Weiss, A. 1992, A&A, 260, 1 

Strigari, L. E., Bullock, J. S., Kaplinghat, M., Diemand, J., Kuhlen, M., & Madau, P. 2007, ApJ, 669, 676, |arXiv:0704.1817| 
Yoo, J., Kochanek, C. S., Falco, E. E., & McLeod, B. A. 2005, ApJ, 626, 51, |arXiv:astro-ph/0502299| 

— . 2006, ApJ, 642, 22, |arXiv:astro-ph/051100l"1 
Zakharov, A. F. 1995, A&A, 293, 1 



